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The LHCb collaboration has studied various promising ways to determine the Unitarity Triangle angle 7. Three 
complementary methods will be considered. The potential of the B — > DK^*> decays has been studied by 
employing the combined Gronau-London-Wyler (GLW) and the Atwood-Dunietz-Soni (ADS) methods, making 
use of a large sample of simulated data. 7 can also be extracted with a time-dependent analysis of Bg ^ DgK 
decays, provided that the Bs mixing phase is measured independently. In addition, the combined measurement 
of the B" — > 7r+7r~ and B^ — > K^K~ time-dependent CP asymmetries allows the determination of 7, up to 
U-spin flavour symmetry breaking corrections. For each method the expected sensitivities to the angle 7 are 
presented. 



1. Introduction 

LHCb aims to study CP violation and rare B-meson 
decays with high precision, using the Large Hadron 
Cohider (LHC) , where aU species of U- mesons are pro- 
duced in 14 TeV pp colhsions [ll, . In these events 
the hh pairs are frequently produced in the same for- 
ward (or backward) direction. The LHCb detector is a 
single-arm spectrometer with a forward coverage from 
10 mrad to 300 mrad in the horizontal plane (i.e., the 
bending plane of the magnet) . The acceptance lies be- 
tween 10-250 mrad in the vertical plane (non-bending 
plane). The detector layout in the bending plane is 
shown in Fig. [TJ 



Muon detector i 




Figure 1: The LHCb setup with the diflterent sub-detectors 
shown in the horizontal (bending) plane. 



2. Event Selection 

LHCb will collect large samples of all types of B 
mesons and baryons. These samples will allow the 
precise measurement of all the three (a, j3 and 7) an- 
gles of the Unitarity triangle and of the Bg mixing 



phase. Here, three complementary methods to extract 
the Unitarity Triangle angle 7 will be considered. 

As it will become clear in the next sections, the 
three complementary methods considered for the ex- 
traction of the Unitarity Triangle angle 7 use differ- 
ent techniques to determine 7. However similar ap- 
proaches are used in the event selection and in the 
selection at the trigger level. In all cases, samples of 
signal and background events are simulated through 
the LHCb apparatus; the selection is then optimized 
in order to maximize the signal efficiency and mini- 
mize the background contribution. The main selec- 
tion criteria used take into account that the large B 
mass produces decay products with high pT and that 
the long B lifetime produces tracks with large impact 
parameter with respect to the primary vertex. Fur- 
thermore particle identification cuts are used for the 
-k/ K identification and cuts on the significance of the 
flight distance of the B are used in order to have a B 
detached from the primary vertex. 

All the following results are for an integrated lumi- 
nosity of 2 fb""'^, which corresponds to one nominal 
year of data taking. 



3. Extracting 7 from B -^ DK decays 

Interfering tree diagrams in the B^ -^ DK^ decays 
allow the determination of the Unitarity Triangle an- 
gle 7. Here, D can be a I?° or a D^ and the D" and 
the Z?" are reconstructed in a common final state. 

In LHCb 3] this is done by employing the combined 
Gronau-London-Wyler (GLW) [^ and the Atwood- 
Dunietz-Soni (ADS) methods. In the first method 
the information of the CP even eigenstates, where a 
D decays to tt+tt" or K~^K~ is used. Note that re- 
construction of the CP odd eigenstates, which include 
neutral particles, is rather challenging in LHCb and is 
not considered in this analysis. The ADS method ex- 
ploits the interference between the favoured and dou- 
ble Cabibbo suppressed decay modes of the neutral D 
mesons to state such as Kit or Ktttttt. Using the ADS 
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method the foUowing equations can be written for the 
decay rates for the neutral B^ -^ D°{^ Kn)K*° 
channels [6| (very similar equation can be written for 
the charged decays, see also [7[) : 

r{B" -^ {K+7r-)DK*") = NKAl + irBTD'f + ^rBrD- 

cos{Sb+Sd+j)) , (1) 

riB° ^ iK-TT+)DK*°) - NiiArl +rl + 2rBrD ■ 

cosiSB-SD+l)) , (2) 

T{B^^{K-^+)dK^) - NKAl + {rBrDf + 2rBrD- 

cos{Sb + Sd- 7)) , (3) 

T{W ^ {K+7:-)dK^) = NkA^I +rl + 2rBrD ■ 

cos{dB -Sd- 7)) , (4) 



3.1. Sensitivity to 7 



where 



As a result of the event selection, the annual yields 
together with the total efficiencies and with the back- 
ground to signal ratios are listed in Table U [g, 0] . 



rB 



\A{B0 ->Wk*o)\ 



is the ratio of the magnitudes between the two ampli- 
tudes of the B-decay and similarly 



Channel 


etot [%] 


S 


Bt-t/S 


B° -^ DO{K+TT-)K*° 


0.33 


3350 


< 2.0 


B°^D°{K+Tr-)K*'' 


0.33 


536 


< 12.8 


B''' ^Dcp{K+K-)K*° 


0.46 


474 


< 4.1 


B° ^ Dcp{n+n-)K''° 


0.36 


134 


< 14 


B+ ^ D'>{K+n-)K+ 


0.50 


28000 


0.63 


B+ ^ DO{K-n+)K+ 


0.50 


100 


7.8 


B+ -^ DO{K+K-)K+ 


0.51 


3000 


1.2 


B+ ^ 00(71+ n-)K+ 


0.58 


1000 


3.6 



Table I Summary of event yield (S), experimental effi- 
ciency (etot), and background to signal ratios. 



ro 



\A{D° -> K+TT- 
\A{W ^ K+n- 



Furthermore Sb represents the strong phase differ- 
ence between the two B-decays while 5d represents 
the strong phase difference between the two D-decays. 
Nkt7 gives the overall normalization and represents 
the total number of B'~' -^ {Ktt)uK*^ events. It can 
be seen that there are two favoured rates (1) and (3) 
and two suppressed rates (2) and (4). Further in- 
formation can be added by including the decays to 
CP-eigenstates, such as tt+tt" or K+K~ and so: 



r{B"^DcpK*") = Ncpil + r% + 
2rBC0s{SB +7)) 



(5) 



r{BO ^ DcpK*°) = Ncpil + rjj + 
2rB cos{Sb - 7)) 



(6) 



Here, Ncp represents the total number of B'^ -^ 
DcpK*^ events. Note that ru is known and is equal 
to 0.060 ± 0.003 [| and Ncp can be calculated from 
Nktt taking into account the different efficiencies and 
branching ratios. Furthermore, for what concerns r^, 
for the neutral B, values smaller than 0.6 with 95% 
probability have been measured at Babar [9|, while 
for the charged B the latest world average value is 
rs ~ 0.10 ± 0.02 |l3]- Sb is unknown for the neu- 
tral i3-meson, while its value is fixed at 130° for the 
charged i?- meson [ll|, [l3 ■ fc is assumed to be in the 
range [-25°; +25°] due to a hmit set by the CLEO-c 
collaboration [13|. In total there are 5 unknowns (7, 
rB, Sb, 5d and Nk-k) and 6 observables. 



Both for the charged and neutral i?-meson decays a 
standalone Monte Carlo (MC) simulation was used to 
generate the event yields and fit the unknown param- 
eters. All the unknown parameters have been scanned 
as shown in Table |TT1 The input value of 7 has been 
fixed at 60°. For the charged i?-meson with one year 



Parameter 


Scan range 




B° ^ D°7^*o 


B± _ D°K+ 


5b 

&D 

rB 


[-180° ; +180°] 

[-25° ; +25°] 

[0.0 ; 0.6] 


fixed at 130° 

[-25° ; +25°] 
[0.0 ; 0.2] 



Table II Scan ranges of the input parameters of the toy 
used for the extraction of 7 for both the charged and the 
neutral _B-mesons. Note that while 5b has been scanned in 
the full range for the neutral B-meson, its value was fixed 
at 130° for the charged B-meson [HI . Il3 ]. 5d has been 
scanned in the range [-25°; +25°] due to a limit set by the 
CLEO-c collaboration [I^. For what concerns rs, for the 
neutral B, values smaller than 0.6 with 95% probability 
have been measured at Babar 9] , while for the charged B 
the latest world average value is rs = 0.10 ± 0.02 [lOj] . 



of data at nominal luminosity (2fb^^), the angle 7 can 
be determined with a precision in the range 8.2°— 9.6°, 
depending on the value of the strong phase 5d [7| • For 
the neutral B the angle 7 can be determined with a 
precision smaller than 10°, depending on the value of 
the strong phase 5b and for r^ values bigger than 0.3 
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4. Extracting 7 from S° -^ DgK decays 

The relations between the B^-meson mass eigen- 
states \Bh,l) and their flavour eigenstates jB"} and 
\B^), can be expressed in terms of linear coefficients 
p and q: 

\BH.L)=p\B'^)TqW) ■ (7) 

The difference in mass and decay rates are defined 
as: 



Af and Aj are the decay amplitudes (e.g. Af 
{f\T\B^)) and 



A/ 



qA^ 

pA, 



qJ2_ 
pAj 



(13) 



For the B^ -^ DfK^ decay channels (see Feynman 
diagrams in Fig. ^ a i?°, as well as a i?°, can decay 
directly to DjK^ or DfK^'. In addition, these rela- 



Am = niH ~ ruL , AF = F^ - F^ . (8) 
The average mass and decay rate are defined as: 

- = — 2 — ' r = ^^— . (9) 

The decay rate at a time t of an originally produced 
B") to a final state / is given by: 



rso^/(t) = l(/|r|B°(t))p 



(10) 



where T is the transition matrix element. The time 
evolution of the flavour eigenstates \B°) and \B^) is 
then given by the four decay equations: 



TB^fit) = \Aff{l + \Xf\^)- 



-rt 



, , AFi , AFt 
(cosh h -D/ smh 

Cf cos Ami — Sf sin Ami) 



B^! 



(t) 



(cosh 



P 

q 
Art 



-rt 



1 + |A/|^ 
+ Df sinh ■ 



2 
AFi 



2 ' 2 

Cf cos Amt + Sf sin Ami) 



.7(^) - 



A 



(l + IAjP) 



-rt 



(cosh 



AFi 



^ . , AFi 

D-F smh 

f 2 



-Cj cos Ami — 5jsin Am,i) 



^B-^jit) = 



^7 



(cosh ■ 



P 

q 

AFi 



(l + IAjP) 



-rt 



2 
AFi 



Dt sinh ■ 
f 2 



-Cj cos Ami -I- 5-Fsin Am,i) 



where 



Df - 
C/ = 
Sf = 



2ReXf 

i-IA/P 

l + IA/P 
2Im\f 



(11) 



(12) 
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Figure 2: Feynman diagrams for the B^ -^ D^ K^ (a), 
bI -^ D+K- (b), B° -^ D+K- (c), bF ^ L>7ii:+ (d) 
decays. 



tions hold for the decay amplitudes: \Af\ = \Aj\ and 
|Aj| = \Af\. Assuming |-| = 1, it is possible to write 
|A/| = |Aj|. The terms A/ and Aj are calculated as 



^D7K+ 



A 



D-K+ 



(v:tyts\ fVubK 



A 



D+K- 



p 



^^^ = lA^- 



D7 K+\ 



yvtbV^ 

,»(A,-(7+0,)) 



v*v 



A 



1^ B, ^DtK- 



Ai 



-lA 



D7 K+\ 



VtbVt*s 

y:byts, 



V^bVcs 



(14) 



where [^2/^1! is the ratio of the hadronic amplitudes, 
which is expected to be of order unity, A^. is the strong 
phase difference between Ai and A2 and 7 + 0s is the 
weak phase. The </><, angle originates from B^ mixing 
and can be measured directly using the B'^ -^ J/'4"P 
decay ^. 

4.1. Sensitivity to 7 + 0^ 

As a result of the event selection, the annual yields 
together with the total efficiencies and with the back- 
ground to signal ratios are listed in Table IIIII 

In a toy MC study multidimensional probability 
functions (PDFs) are constructed. These are sup- 
posed to mimic the outcome of an analysis of data 
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S 


etot [%] 


Bti/S 


B° -^ DfK^ 


6200 


0.32 


< 0.7 



Table III Summary of event yield (S), experimental effi- 
ciency i^tot), and background to signal ratios. 





Anisips-^) 


|A/I 


^Tl/T2 (°) 


7 + 0. (°) 


Input value 


17.5 


0.37 





60 


Fitted value 


17.5 


0.372 





60.4 


^(ly) 


0.007 


0.061 


10.3 


10.3 


'^(5y) 


0.003 


0.027 


4.6 


4.6 



Table V Summary of the sensitivity results for the toy for 
B° -> DfK"^ and B° ^ D7-K+ events. 



acquired at LHCb. So first the behaviour of the ex- 
periment is described building different PDFs. In the 
toy in use for the sensitivity studies of the B'^ ^ DgK 
channels several PDFs describing the mass distribu- 
tion, the proper time acceptance, the flavour tag- 
ging and the particle identification response have been 
built both for the signal events and for the background 
events according to the studies done with the full 
Geant 4 simulation [15]. The final PDF is built as 
the product of all the different PDFs and finally a 
likelihood fit to the generated data is performed to 
extract the parameters and their errors. In this toy 
not only the B^ — > DgK decay channels are consid- 
ered, but also the topologically similar S" -^ DsTt 
decay channels. The B'^ ^ Dgir decay channel has a 
larger branching ratio and consequently a higher an- 
nual yield (140000 events) and it allows for the deter- 
mination of the Arris parameter. A summary of the 
input parameters that are used is given in Tab. IIVI 
The final sensitivity results are shown in Tab. |Vl It 



5. Extracting 7 from 5° 

B2 -^ K+K- 



n ' TT 



and 



In this section the way to extract 7 through the 
combined measurement of the B'^ — > 7r"'"7r~ and B^ -^ 
K^K~ CP asymmetries and under the assumption of 
invariance of the strong interaction under the d and s 
quarks exchange (U-spin symmetry) (l7| is described. 
In Fig. [3] the SP x — > h'^h'~ tree diagrams are shown. 



d{s) 



fDO\ 



B^ B 




^- {K- 



dis) 



d{s) 



Figure 3: The B,^-. -^ h h' tree diagram. 



Parameter 


Input value 


cr(mso) (MeV) 


14 


Ar,/r, 


0.1 


Ams{ps~'^) 


17.5 


mistag fraction lj 


0.328 


tagging efficiency etag 


0.5812 


|A/| 


0.37 


7 + 0. (°) 


60 


Ati/T2 (°) 





B^ — > D~n~^ annual yield 


140k 


Bs -^ DjK^ annual yield 


6.2k 


B° -^ DJ-K+ B/S ratio 


0.2 


B° -^ DfK^ B/S ratio 


0.7 



Table IV Input parameter values for the toy MC simula- 
tion program. 



can be seen that with one year of data taking at nom- 
inal luminosity a sensitivity on 7 -I- ^^ of 10.3° can be 
obtained. A more detailed study can be fomid in fig. 



For a neutral i3-meson decaying into a CP eigen- 
state /, the time-dependent CP asymmetry is given 
by: 



Acp{t) 



nB%/s{t) 
nB^d/sit) 



■f)-nB',/sit)^f) 



+ 

(15) 



—Ccp cos Ami + Sep sin Amt 
" cosh ^t-A^lsm\i^ 



'-t 



where r(B"<i/,(i) ^ /) and r(B;]/^(t) ^ /) are the 

decay rates of the initial B and B states respectively, 
and Am and AP are the mass and width differences 
between the two mass eigenstates. 

As shown in [13] Ccp and Sep for the B^ -^ tt+tt^ 
and B'l -^ K^K~ can be written as functions of the 
Unitarity Triangle angle 7, and of the two hadronic 
parameters d and (which parametrize respectively 
the magnitude and phase of the penguin-to-tree am- 
plitude ratio), as follows 






2(isin0sin7 



\ — 2d cos 6 cos 7 + d^ 

sin {(l)d + 27) — 2d cos 9 sin {(pd -I- 7) -I- d^ sin cjid 



1 — 2d cos 6 cos 7 + d^ 
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Ckk 
Skk 



sin 



2d' suvO' sin 7 
1 + 2d' cos 0' cos 7 - 

(0s + 7) + d'"^ sin ( 

~dF^ 



-d'2 



[(f)s + 27) + 2d' cos 9' sin (0^ 
1 + 2d' cos 6*' cos 7 - 



where (j)d and (/)s are the B'^/B'^ and B^/B^ mix- 
ing phases which in LHCb wiU be measured from 
i?° — > JN)Ks decay and from i3° ^ ^/V"/* respec- 
tively |14l |. This system of four equations and five 
unknowns {d,d',9,0' and 7) can be solved with the 
help of the U-spin symmetry, as a consequence d = d' 
and 9^9'. This results in an over-constrained system 
of three unknowns and four equations. 



5.1. Sensitivity to the CP asymmetries 
and to 7 

As a result of the event selection, the annual yields 
together with the total efficiencies and with the back- 
ground to signal ratios are listed in Table IVTl 



Event type 


BR 

(xlO-«) 


S 


etot [%] 


Bti/S 


Bsp/S 


BS - n+n- 


4.8 


35700 


0.93 


0.46 


0.08 


Bl - K^.- 


18.5 


137600 


0.93 


0.14 


0.02 


B'l -. n+K- 


4.8 


9800 


1.02 


1.92 


0.54 


B'' ^ K+K- 


18.5 


35900 


0.97 


< 0.06 


0.06 



Table VI Summary of event yields (S), Branching ratios, 
experimental efficiency {etot), and background to signal 
ratios, considering both the bb combinatorial background 
and the specific background. 



As for the B^ ^ DgK decays a toy MC is used to 
mimic the outcome of an analysis of data acquired at 
LHCb. In this toy, used for the CP sensitivity studies 
of the B9. — > h~^h'~ channels, several PDFs describ- 
ing the mass distribution, the proper time acceptance, 
the flavour tagging and the particle ID response have 
been built both for the signal events and for the back- 
ground events according to the studies done with the 
fuU Geant 4 simulation [18^ . The final PDF is buih as 
the product of all the different PDFs and a likelihood 
fit to the generated data is performed to extract the 
parameters and their errors. 

The extraction of 7 is then performed using a 
Bayesian approach in three different U-spin scenarios: 



• Assuming perfect U-spin symmetry: d ~ d' and 
9 = 9' . 

• With a weaker assumption on the U-spin sym- 
metry: d = d' and no constraint on 9 and 9' . 



• With an even weaker assumption on the U-spin 
symmetry: ^ = d'/d = [0.8,1.2] and no con- 
straint on 9 and 0' . 



In the three different cases the sensitivity on 7 is taken 
^ -'^ainsidering the 68% probability interval of the result- 
ing PDF distribution for 7. 

Fig. |4] shows an example of a resulting PDF for 7 
obtained assuming perfect U-spin symmetry. The 68% 
and 95% probability intervals are visible. The 68% 
probability interval corresponds to a sensitivity of 4° . 
The sensitivity in the second and in the third scenarios 
increases and varies between 7° and 10°. 



>, 0.015 




y[°] 



Figure 4: 7 probability density obtained in case of perfect 
U-spin symmetry. The 68% and 95% probability intervals 
are visible. The 68% probability interval corresponds to a 
sensitivity of 4°. The plot has been obtained using tools 
developed by the UTfit Collaboration [l9l ]. 



It is important to note that the extraction of 7 by 
means of the i?° — > tt+tt" and B^ -^ K^K~ decays 
uses not only tree diagrams but also loop diagrams 
and is therefore sensitive to new physics. 



6. Conclusions 

In these proceedings three complementary methods 
for the extraction of the Unitarity Triangle angle 7 
have been discussed. 

The potential of the B -^ DK'^*^ decays has been 
studied by employing the combined Gronau-London- 
Wyler (GLW) and the Atwood-Dunictz-Soni (ADS) 
methods. For the charged B-meson, with one year of 
data at nominal luminosity, the angle 7 can be de- 
termined with a precision in the range 8° — 10°, de- 
pending on the value of the strong phase 6d 01 ■ For 
the neutral B, the angle 7 can be determined with a 
precision smaller than 10°, depending on the value of 
the strong phase 5b and for rs values bigger than 0.3 
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It has been shown that the angle 7 can also be ex- 
tracted with a time-dependent analysis through the 
Bg — » DgK decays, provided that the B^ mixing phase 
is measured independently. With one year of data tak- 
ing at nominal luminosity, a sensitivity on 7 + 0s of 
10° can be obtained. 

The combined measurement of the B^ — > ■k~^'k~ and 
B^ -^ K^K~ time-dependent CP asymmetries allows 
the determination of the Unitarity Triangle angle 7, 
up to U-spin flavour symmetry breaking corrections. 
Here a sensitivity of 10°, with one year of data taking 
at nominal luminosity, can also be obtained, but the 
final results depend on the assumption on the breaking 
of the U-spin symmetry. This method uses not only 
tree diagrams but also loop diagrams and is therefore 
sensitive to new physics. 

Other methods, including B -^ D{KKTm)K and 
B -^ D{KaTT~^Tr~)K decays, to extract 7 at LHCb 
have been studied. More details can be found in '20'. 

El. 
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